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ABSTRACT 

Context. The sterile neutrino is an excellent dark matter candidate, which can be searched for in a wide range of astrophysical sites. 
It has previously been shown that the optimal search strategy is to consider dwarf galaxies belonging to the Milky Way. 
Aims. We search for line emission from decaying dark matter. 

Methods. We analyse publicly available Chandra X-ray observations of the dwarf galaxy Draco. 

Results. The Draco and blank sky (background) spectra are nearly identical in shape, which allows us to conclude: i) Dwarf 
spheroidals are ideal for studying dark matter X-ray emission since the baryonic noise is impressively low, ii) there is very little 
room for line emission, which leads to constraints in the mass-mixing angle parameter space of the sterile neutrino. We compare the 
standard flux derivation method to a very conservative rebinning approach. The resulting constraints are strongly dependent on the 
chosen method. 
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1. Introduction 

The sterile neutrino is a strong particle candidate for dark 
matter (DM). With just three sterile neutrinos (gauge sin- 
glets), one easily obtains the correct abundance of dark 
matter, a very simple explanation for the observed fla- 
vor oscillations of the active neutrinos as well as their 
mass splitting, and, in additio n, is a natural explana- 
tion for the baryon asymmetry ( Dodelson & Widrowl IT 994t 



Shi & Fulled 1 1 999h iDolgoy & Hansen! 120021: lAbazaiian et al 



2001 at lAsaka et alJ 120051: lAsaka & Shaposhnikovl 120051) . 
The underlying particle model, called t he vMSM, is de- 
scribed in detail in a number of papers (lAsaka et all 120051: 
Asaka & Shaposhnikov 2005; Asaka et al. 2007; Shaposhnikov 
2007; Gorbunov & Shaposhnikov 2007; Laine & Shaposhnikov 
2008: lShaposhnikovl l2008h. and also in the excellent re- 
view by iBovarskv et al.l d2009bl) . Additionally, the sterile 
neutrino may have interesting effects on a range of differ- 
ent astrophysical objects, including as an explanation for 
pulsar kick velocities, facilitating core collapse supernova 
explosions, affecting early star formation, reionization and 
struct ur e formation, or a ss isting in flation dKusenko & S egre 
19971: iFrver & Kusenkol 120061: lHidaka& Fulled |2006 
Biermann & Kusenko 2006; Hansen & Haiman 
Mapelli et al. 2006; Bezrukov & Shaposhnikov 
Shaposhni kov & Tkachevj 120061; [K usenko et al.l 
Petraki & Kusenkol l2008t iPetrakii 120081: iBovanovskyi 
Gorbunov et alJ 120081) (see IBovarskv et alJ d2009b) ' 



2004 
2008 
2008; 
2008; 
for an 
extensive list of references). 

The large range of astrophysical objects affected by the ster- 
ile neutrino allows for independent measurements of these new 
particles. The two free parameters for the lightest sterile neu- 
trino, the DM particle, is its mass, m s , and mixing angle with the 



active neutrinos, sin~(20), and various observations have already 
excluded large parts of this parameter space. 

Structure formation will be changed if the DM particle 
is warm as compared to cold, and for instance Lyman-a ab- 
sorption in quasar spectr a allows one to place lower bounds 
on th e DM particle mass (lHansen et all 120021: IVM et all 120061, 
2005 : ISeliak et al] 12006 ). The most recent Lyman-a- analysis 
(Bov arskv et al. I l2008bld) finds that for all masses above m s = 
2 keV, there exist parts of the parameter space in agreement with 
all observations. The originally p roposed particle productio n 
method, non-resonant production (Dodelson & Widrow 1994), 
is already rule d out observationall y, and only the resonant pro- 
duction (RP) dShi & Fulled fl999). with a large initial lepton 
asymmetry, is in agreement with all observations. This lepton 
asymmetry can be constrained directly fr om nucleosynthesis 
(ISerpico & RaffeU2005l:lDolgov et al.l2002l) . Given this limit on 
the lepton asymmetry the RP needs a sufficiently large mixing 
angle to make the abundance of the sterile neutrino as large as 
the observed DM abundance. This gives a strong lower limit on 
the mixing angle, which is approximately sin 2 (20) > 10" 12 for 
m s = 2 keV dLaine & Shaposhnikovl2008l) . 

A different, but very firm, lower boundary on the mass 
is obtained through the phase space density of nearby dwarf 
galaxies. The Tremaine-Gunn bound dTremaine & Gunnlll979l) 
gives a model indep endent boundary of roughly 0.4 keV 
(Bovarskv et al. 2009a). This boundry can be strengthened if the 
production method is known, and for the R P the boundary is ap- 
proximately 1 keV (IBovarskv et al.ll2009al) . 

The first constraints on the mixing angle came fr om the non- 
observation of decay lines (Dolgov & Hansen 2002). This is be- 
cause the sterile neutrino decays through a 1 -loop into a photon 
and an active neutrino, w ith a decay rate of dPal & Wolfensteinl 
H98l iBarger et1d1l995h : 
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Since the DM decay almost at rest, the decay line is very 
narrow and e asily searched for in X-r a y and soft gamma ray 
observations (Dolgov & Hansen 2002; Riemer-S0rensen et al. 
2007t lAbazaiian et all l 2007t faemer-Sdrensen et alJ 12006: 



Bovarskv et al. 



Bovarskv et al. 



2006b, 
2006a, 



2007b; Abazaiia netai] 



200831 lYiiksel et al l 



Wats on et al.ll2006l:lLoewenstein et alJ fcoOS ) 




The main problem with X-ray observations is the bary- 
onic background, and one therefore needs to examine sites 
with very few baryons. Such systems alread y considered 
include the dark blobs in merging clusters (lHansen et al.l 
l2002t lRiemer-S0rensen et alJ 120071 lAbazajian et all 120071) or 
simply looking out through our own Milky Way halo 
dRiemer-S0rensen et al.ll200"6tlBovarskv et a l. 2006b). However, 
it was demonstrated that the optimal place to search for the de- 
cay line of the sterile neutrino is from the dw arf galaxies of 
our Milky Way such as Ursa Minor, Draco etc. (iBovarskv et al.l 
2006b, 2007b). These dwarfs are nearby, they are X-ray quiet, 
have a high central DM density, and at the same time we 
have fairly accurate mass models from optical observations 
dGilmore et al.ll2007t IWilkinson et"aT1l2004t IStrigari et ai1 l2008 , 
120071) . 

In this paper, we analyse Chandra X-ray data of Draco. The 
Draco and blank sky spectra are almost identical and leave no 
significant room for line emission from decaying DM. The ab- 
sence of any line signal is used to derive constraints in the mass- 
mixing angle parameter space of the sterile neutrino. 



2. The observations 

2. 1 . X-ray analysis 

There exist two public Chandra observations of Draco with ob- 
servation ids 9568 (24.8 ks) and 9776 (12.2 ks) in the NASA 
HEASARC archival Before extracting spectra, we reprocessed 
the data with the newest versions of CIAO (4.1) and CA LDB 
(4.1) following standard procedures dFruscione e t al. 2006)). A 
total of 12 evenly distributed point sources were removed and 
the light curve cleaned, reducing the total exposures to 20.3 ks 
(9568) and 12.0 ks (9776)). The central part of Draco is situ- 
ated on the S3 chip. The spectra were extracted from a square 
region of (7.6 ') 2 ~ (0.18kpc) 2 , avoiding the edges of the chip. 
The response matrices were generated and the spectra and re- 
sponse matrices of the two observations were combined using 
/foo/J3 (justified by the very similar count rates in the two oth- 
erwise identical observations). For background subtraction we 
extracted the spectrum of the identical region of the chip from 
the blank sky data provided with CIAO. 

The Draco and blank sky spectra are shown in Figure [1] 
where it is clearly seen that they are almost identical, with very 
little room for a dark matter signal. 

The normalization of the blank sky spectrum was determined 
from the average count rate in the 4-6 keV interval for both 
blank sky and observation. This interval was chosen because 
both spectra are very flat and well calibrated here. Additionally, 
no thermal gas emission is expected at these energies, as Draco 
is a small and cold structure. The uncertainty of the average 
count rate is less than 0.5%. Normalising at higher energies gives 
slightly different background levels e.g. 10-12 keV gives a 3.5% 
lower level, and 12-14 keV (which is usually used for galaxy 
cluster observations) gives a 0.5% higher level. However, both 
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Fig. 1. The spectrum of the central (7.6 ') 2 (black) and the corre- 
sponding normalized blank sky spectrum (green). 
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1 http://heasarc.gsfc.nasa.gov/docs/archive.html 

2 ftools, HeaSoft, http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftoolsaMteiilchlM a monochromatic emission line 
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Fig. 2. The background subtracted spectra of the S3 (center of 
Draco, fat black), S2 (8.6 ' from center, green), and S3 (17.3 ' 
from center blue) chips rebinned for visualisation. In the case of 
any observable signal from dark matter, there will be a radial de- 
pendence, like the one around 1 keV, where the signal disappears 
far from the center. 



intervals have higher uncertainties on the normalisation due to 
the large scatter. 

The background-subtracted Draco spectrum is shown in 
Figure|2](black). Apart from a small signal »1 keV, the spectrum 
is effectively zero. From this we see, i) that dwarf spheroidals 
have no significant baryonic X-ray emission, ii) there is virtu- 
ally no room for a dark matter decay line. 

After the background subtraction, the only possible signal is 
at low energies (« 0.7 - 1.5 keV). The best fitting Gaussian has 
a width of > 0.5 keV, which is an order of magnitude larger than 
the width of the instrumental resolution at 1 keV. To check for 
any radial dependence in the signal, we extracted spectra from 
similar regions of the S2 and SI chips centered « 8.6 ' and 
« 17.3 ' from the center of Draco. The three spectra are com- 
pared in Figure |2] The low energy "bump" around 1 keV is not 
clearly present in the outer radial bins. It is tempting to believe 
that this is a signal proportional to the mass within the field of 
view. Nonetheless, the excess is not very significant and is too 
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Fig. 3. The rebinning approach for E y = 5 keV. All data within Fig. 4. The Draco spectrum and the background model (green 
the interval [E y - E FWHM /2 : E y + E FWHM /2] was rebinned to line) including the potential line at E y = 5 keV. For visualisation, 



one bin (central data point). The width of the Gaussian is given 
by the resolution and the maximum by the 3<x upper limit of the 
rebinned data. 



backgrounds such as stellar and thermal radiation will also fol- 
low the mass distribution of Draco. 

The upper limit on the flux was derived in the 0.7-10 keV 
interval using the sp ectral fitting package Xspec (version 12.4, 
lArnaud et alj 12007) following different approaches. The flux 
derivations are based on the dark matter in Draco being practi- 
cally at rest (v/c ss 10~ 4 . IGilmore et af]|2007D so the line broad- 
ening due to internal motion is negligible compared to the in- 
strumental resolution. The instru mental resolution is given by: 
Efwh m - O.Q\2E y + 0.12keV dProposers Observatory Guide! 
l2007h . 

One employed approach is to rebin the data with a binwidth 
corresponding to the instrumental resolution. At a given energy, 
E y , all bins within the interval [E y - E F whm /2 : E y + E F whm /2] 
were rebinned to give a single bin value and uncertainty (fat data 
point in Figure [3j. The possible line emission was defined as 
a Gaussian with the width given by the resolution (EfWHM = 
2.35£ , cr for a Gaussian distribution) and maximum at the 3cr up- 
per limit of the rebinned data. The upper limit on the flux was 
determined from the Gaussian over the same interval as the data 
were rebinned. 

Another approach was to model the back ground instead 
of sub tracting it (same method as discussed in Bovarsk v et alj 
2006b). The background was fitted with a continuum consisting 
of an exponential and 10 Gaussians. The fit gave a reduced^ 2 of 
1.3 (for 637 d.o.f.). The line emission was added to the model as 
a Gaussian with the width given by the instrumental resolution 
(see Figure @}. The total model was fitted to the data with the 
normalisation of the Gaussian as the only free parameter. The 
flux was determined exclusively from the Gaussian with the nor- 
malisation set to the 3cr upper limit from the fit. 



2.2. Observed mass 

The dark matter mass within the field of view is determined from 
th e density profile deriv ed from photometric and kinematic data 
by IGilmore et al.l(l2007l) . The density profile in the 0.1 - 0.5 kpc 
range is well described by p(r) - po(l + (r/ro) 1 ' 5 ) -2 ' 25 where 
po = 0.65 x 1O 9 M kpc 3 and r = 0.28 kpc. We extrapolated 
this profile to the observed range of - 0.18 kpc and integrated 



the line normalisation has been increased by a factor of 10. 
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Fig. 5. The obtained upper limits on the line emission flux. 
The modelling Gaussian approach (solid red) gives stronger 
constraints than the conservative rebinning approach (black 
hatches). 



along the line of sight to determine the mass within the field 
of view. The t otal dark matter mass of Draco within 0.5 kpc is 
> 6 x 10 7 M Q (iGilmore et al.ll2007HStrigari et al.ll2008l) of which 
« 6.7 x 10 6 M Q is within the observed square. The statistical un- 
certainty on the ma ss is very small, but the systematics are ver y 
hard to determine (IGilmore et al.ll2007t iBovarskv et al.ll2007bl) . 



3. Results 

Figure|5]shows the obtained upper limits on line emission flux as 
a function of photon energy (coloured areas are excluded). These 
constraints are very general and apply to all dark matter candi- 
dates with a monochromatic line emission. The vertical lines in 
the rebinned approach constraints are single energies, where the 
upper limit on the flux becomes exactly zero for a 3cr upper limit 
on the flux (using this approach with the chosen background nor- 
malisation). If we require e.g. 5cr, they will disappear. 

The fluxes are converted to constraints in the m s - sin 2 (20) 
parameter space for sterile neutrinos of the Majorana type, as- 
suming the sterile neutrinos to account for all of the dark matter 
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Fig. 6. The obtained constraints in the m s — sin 2 (26) parameter 
space for sterile neutrinos of the Majorana type. The modelling 
Gaussian approach (solid red) gives stronger constraints than 
the conservative rebinn ing approach (black h atched). Earlier 
constraints (taken from Bovarsk v et alj l2009bl and references 
therein): Phase space considerations (grey below 1 keV), nucle- 
osynthesis (grey region below dotted line), non-resonant produc- 
tion giving Q. > Q.DM (grey region above dashed line), earlier 
observational constraints (blue). 



in Draco (lRiemer-S0rensen et al.lf2 006; Bo varskv et al 1200 6b): 



sin 2 (2<9) < 1 x 10 1 



Fdet 



ergcm 2 s 1 



IkeV/ 



(M fov /M e ) 



(D L /Mpc) 2 



-l 
(2) 



The distance to Draco is 80 kpc (iGilmore et alJl2007h . The re- 
sulting conservative constraints are shown in Figure [6] (red and 
black hatched) compared to earlier constraints fr om galaxy clus- 
ters, g alaxies, and the Milky Way halo (blue, iBoyarsky et al.l 
2009b, and references therein). The overlap with earlier con- 
straints is an excellent confirmation of the excluded parameter 
space, since the Draco constraints are complementary and in- 
dependent of earlier constraints. Even with the relatively short 
exposure, the constraints are good due to the almost perfect sim- 
ilarity of the spectra of Draco and the background. The differ- 
ence between the Gaussian modelling (red) and the rebinning 
(black hatched) approaches is clearly visible, with the Gaussian 
modelling method providing stronger constraints at most ener- 
gies than the more conservative rebinning method. 



4. Future observations? 

The limiting factor of the constraints are the uncertainties in the 
observations and blank sky, which can be improved by increas- 
ing the exposure time. The statistics improve with the square 
root of the exposure time so a 100 ks observation improves the 
constraints by a factor of « 2. In order to exclude the entire cur- 
rently allowed region, the needed exposure time for Chandra 
or similar instruments becomes unreas onably high and a dif - 
ferent instrumental approach is needed (Bovarskv et al. 2007a). 
One option could be gratings, but unfortunately gratings are 
not useful for extended sources because of smearing due to 
the ambiguity between energy and incident angle of the photon 
dRiemer-S0rensen et a l. 2007). 



5. Summary 

We have searched for line emission from decaying dark mat- 
ter particles in Chandra X-ray observations of the Milky Way 
dwarf galaxy Draco. The Draco and blank sky background spec- 
tra have a nearly identical shape. This confirms dwarfs as ideal 
for studying dark matter emission since the baryonic contamina- 
tion is close to zero. No obvious line signal is found, which leads 
to conservative constraints in the m s - sin 2 (2#) parameter space 
for sterile neutrinos. Longer exposures are needed to reduce the 
statistical uncertainty and improve the constraints. 
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